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Abstract. - 

We consider three spinless fermions free to move on 2d square lattice with periodic boundary 
conditions and interacting via a U/r Coulomb repulsion. When the Coulomb energy to kinetic 
energy ratio r s is large, a rigid Wigner molecule is formed. As r„ decreases, we show that 
melting proceeds via an intermediate regime where a floppy two particle molecule coexists with 
a partially delocalized particle. A simple ansatz is given to describe the ground state of this 
mesoscopic solid-liquid regime. 



Ordered arrays of charged particles with long range Coulomb repulsion have been a continu- 
ous subject of interest in various branches of physics, including colloidal suspensions, ion rings, 
atomic-ion Wigner crystals, quantum computers, biophysics, plasmas, electrons deposited on 
liquid Helium surfaces, charges created in semiconductor or organic field effect devices. These 
arrays can melt, exhibiting a transition from collective towards independent-particle motion, 
either as a function of the temperature (classical melting) or as a function of the charge 
density (quantum melting) at very low temperature. In principle, the quantum melting can 
be observed using electrons trapped in quantum dots || |J or cooled ions confined in radio 
frequency traps H. Very often, a parabolic confinement is imposed. When the confinement 
is weak and at a sufficiently low temperature, the Coulomb repulsion dominates the kinetic 
energy, the charges are ordered and a Wigner molecule is formed. If the confinement becomes 
stronger, the kinetic energy dominates the Coulomb repulsion, the molecule melts and one 
gets a Fermi system of weakly interacting particles. In a parabolic 2d trap, the molecule 
consists of well-separated shells. Both for the classical melting fl (achieved by increasing the 
temperature for a given trap) and for the quantum melting Q (achieved at zero temperature 
by reducing the size of the trap), it has been observed that melting proceeds in two stages: 
first neighboring shells may rotate relative to each other while retaining their internal order, 
second the shell broadening leads to radial melting. Wigner quantum crystallization in 2d 
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electron dots is characterized by two distinct - radial and orientational - ordering transitions. 
However, a parabolic trap does not yield a uniform charge density, the low density shells at 
the edges could order before the high density part in the bulk, and this two stage melting 
could be an artifact due to the interplay between surface and bulk orderings. It is therefore 
interesting to study if a multi stage melting persists in a system of uniform charge density, for 
instance when the charges are confined on a 2d torus. One has then to take into account the 
translational symmetry of a 2d torus instead of the rotational symmetry of a parabolic trap. 

Another important issue can be mentioned, assuming that insights gained through inves- 
tigations limited to small systems provide the foundations for understanding larger systems. 
Long ago, it was suggested by Andreev and Lifshitz [|| that low temperature localized defects 
change into excitations that move practically freely through a crystal. As a result, the number 
of sites of a quantum crystal may be smaller than the total number of particles present in 
the system for intermediate couplings, such a crystal being neither a solid, nor a liquid. Two 
kinds of motion are possible in it; one possesses the properties of motion in an elastic solid, 
the second possesses the properties of motion in a liquid. 

An intermediate regime of melting was detected using N — 4 spinless fermions in a L x L 
lattice, and it was observed |M |) that a combination of a few plane waves and site orbitals 
was able to describe it, suggesting a liquid-solid regime consistent either with a scenario a la 
Andreev-Lifshitz, or with a possible quantum liquid crystal regime || [To| . This intermediate 
regime was shifted to lower ratios r s when site disorder was included. Many signatures of a 
novel ground state (GS) were observed for intermediate couplings, considering the structure 
[ p"l| of the persistent currents when the torus was pierced by an Aharonov-Bohm flux, the 
statistics J9| of its low energy excitations, the failure B of the Hartree Fock approximation to 
describe the persistent currents above a first ratio r^, the suppression of the same currents 
above a higher ratio rY , the GS magnetization |T^] when the spin degrees of freedom where 
included. Moreover, it was noticed in Ref. |llj that the ratios rf < r s < where the novel 
mesoscopic GS was observed were consistent with those where transport measurements using 
dilute electron gases in 2d field effect devices ]l3[ indicate the puzzling possibility of a novel 
2d metal. 

The purpose of this work is to reveal the exact nature of the intermediate mesoscopic GS 
and to describe it with a simple wave function, considering N = 3 spinless fermions with U/r 
Coulomb repulsion in a L x L square lattice with periodic boundary conditions. Using the 
operators cj (cj) , Ck (c k ) which create (destroy) a spinless fermion either at the lattice site 
j = (jxijy) or in a plane wave state of momentum k = (k x , k v ) of this lattice, the Hamiltonian 
reads: 

H = "* E C h' + \ E 7?T = E e *4 c fc + E U(q)cl +q cl_ q c k ,c k . (1) 

j<j' M k k,k' ,q 

rij — CjCj, dj ji is the shortest distance between sites j and j', — —2t(cosk x + cosk y ) 
and U(q) = U/(2L 2 )J2-cos(qj)/d ^j- The Coulomb energy to kinetic energy ratio r s = 
U j (2ty/Trn e ) for a filling factor n e = N/ L 2 . The operators Ck and Cj are related by the usual 
Fourier transform: 

n yX< (2) 

3 

In the eigenbasis of the non interacting system (eigenvectors cj^c^c^ |0), |0) being the 
vacuum state), the Hamiltonian matrix is block diagonal, each block being characterized by 
the same conserved total momentum K = k%. Moreover, only the non interacting states 
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Fig. 1. - Left sides: Scheme of two low energy Wigner molecules in the limit r s — > oo and of total 
momenta K for L — 8. Right sides: Corresponding relative fluctuations Ad/ {d) of the three particle 
spacings as a function of r s . Smallest (continuous line), intermediate (dashed dot line) and largest 
(dashed line) spacings. The thin dotted lines give the t/U perturbative behaviors. The arrows indicate 
the two crossover ratios r* 



having in common one k out of three can be coupled by the interaction inside a K sub-block. 
Therefore, each K sub-block is a sparse matrix which can be exactly diagonalized using the 
Lanczos algorithm when L is small enough. 

In the large coupling limit (r s — > oo) the eigenstates correspond to rigid Wigner molecules. 
For a localized center of mass, the charges would be simply located on three lattice sites a, b 
and c, the location of those sites minimizing the electrostatic energy. However translational 
invariance implies a derealization of the center of mass in a plane wave state of momentum 
K, and the Wigner molecule wave functions become: 

|*) = i^e^) C t +aC t +6 ct +c | ). (3) 
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For a given shape (a, b, c) of the three particle molecule, one has three well defined inter particle 
spacings d min < d int < d max . The electrostatic energy becomes U{d^ ln +d^ t +d m 1 ax ) while the 
kinetic energy cx t e ff (center of mass effective band width) — > as r s — > oo. For some arbitrary 
values of K and for some molecular shapes of low electrostatic energy, we have decreased r s 
and followed the corresponding levels, ignoring possible level crossings with other levels of 
different K and of different molecular shape. Two examples calculated using a L — 8 lattice 
are shown in Fig. 1, revealing the generic scenario for the melting of a three particle Wigner 
molecule. If one considers the change of the relative fluctuations Ad/ (d) as one increases r s , 
one can see both for d m i n and d max a crossover from a weak coupling behavior where the 
fluctuations are large (the d are not well defined) towards a large coupling behavior where the 
fluctuations become negligible (the d become well defined) . The weak (strong) coupling limits 
can be described using U/t (t/U) perturbative expansions. For instance, the large coupling 
behavior of the Wigner molecule is given at first order of a t/U expansion by: 

l*(l)> = l*> + E^I*a>, 

where the |\I/ Q ) label the Wigner molecules of same K obtained by moving one of the sites 
a, b, c of |\&) by one lattice spacing, AE a cx U being the corresponding changes of electrostatic 
energy. This gives the t/U decays of the three Ad/ (d) indicated in Fig. 1. 

The main point to notice is the clear separation between the crossover ratios r* (indicated 
by the arrows in Fig. 1) characterizing d m i n and d max . As one can see in the data, there are 
relatively large intervals of intermediate couplings where d m i n is well defined while d max is 
not, giving rise to an intermediate behavior for di n t- This tells us that the generic melting of 
a three particle molecule proceeds also in two stages, if one considers a 2d system of uniform 
density, as it had been shown using 2d parabolic traps. The intermediate regime of melting 
consists of a floppy two particle molecule co-existing with a third delocalized particle. 

We now study the ground states (GSs) of the three body problem. The GS momenta and 
degeneracies depend on L, as the possible existence of GS level crossings. For simplicity, let 
us consider the case where L is even. At U = 0, one has a sixfold GS degeneracy which 
is partially removed by an infinitesimal U, the energy of a set of four states with momenta 
K = (±2n/L, ±2tt/L) and K = (±2n/L, T^/L) becoming different to those of the two K = 
states. Using a U/t expansion, one finds that the GSs remain in the first set for L < 6 while 
they go into the second set for L > 8. At t = 0, the low energy Wigner molecules are L? 
triangles (d m i n — dmt = L/2,d max = L/y/2) having L 2 /A different locations of their centers 
of mass and 4 different orientations. This L 2 degeneracy is removed when one turns on t. The 
energies Eq{K) of the L 2 first levels are given when t/U — > by: 

E (K) fa A — 2t eff (cosK x + cosK y ) + 2r eff (cos(K x L/2) + cos(K y L/2)) , (5) 

where A is a K-independent energy, t e ff <x t(t/t r ) JV_1 is the effective center of mass band 
width while r e // cx t{t/U) L l' 1 ^ 1 is the effective band width coming from single particle motions 
which couple triangles of same center of mass but of different orientations (L/2 one particle 
hops). For L > 8, t e fj » r e ff, one has a non degenerate K = GS when t/U — > oo, and 
consequently a GS level crossing as r s increases inside the K — subspace between the two 
weak coupling GSs and the single large coupling GS. For L = Q,r e ff and t e f / are both cx t 3 /U 2 , 
the GSs keep as r s varies their momenta K = (±2n/L, ±2tt/L) and K — (±2tt / L , ^p2ir / L) 
(fourfold degeneracy) and no GS level crossing occurs. Hereafter we study the L = 6 GS of 
momentum K = (2tt/Q,2tt/6). This allows us to avoid the complications coming from the 
GS level crossing for L > 8. However, as shown by the previous examples, our results will 
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Fig. 2. - Left figure: Scheme of a x-oriented partially melted triangular molecule (x-PMTM) for 
L — 6. The arrows give configurations allowed by a partial t/U expansion of the 2PWM. Right figure: 
GS projections P x (K,k x — 0) (solid line) P x (K,k x = 2n/3) (dashed-dotted line), P x (K,k x = — 2ir/3) 
(dotted line) over x-PMTMs of momentum k x . P(K, k x = 0, k y = 0) (upper thin dashed line) gives the 
total projection over the PMTM of momenta k x = and k y = 0. Total momentum K = (27r/6, 2n /6). 



be relevant to generically describe the multi stage quantum melting of a N = 3 low energy 
Wigner molecule if one continuously follows a given level from large couplings towards weak 
couplings. 

For L = 6, the degeneracy of the L 2 triangular molecules is broken when one turns on t. A 
t/U expansion gives for the L 2 = 36 low energy levels 

t 4 

Eq{K) =A 2 - 2< 3 (cos 2irK x + cos 2nK y ) + 2r 3 (cos 2n3K x + cos 2tt3^) + 0(^ ).. (6) 

A 2 = 0.9023L/- 208.9(t 2 /U), t 3 = 1000{t 3 /U 2 ) and r 3 = I660(t 3 /U 2 ). This t/U expansion 
makes sense when both d m i n and d max are well defined, with small relative fluctuations of 
order t/U. This means r s > 200 (see Fig. ||). To describe lower r s (40 < r s < 200), where 
dmin ~ 3 is well defined, while d max has still large fluctuations, we propose a simple ansatz 
based on the concept of partially melted triangular molecules (PMTMs) . A x-oriented PMTM 
(x-PMTM) is a rigid two particle Wigner molecule (2PWM) with d m in — L/2 combined with 
a third particle free to move with a wave vector k x parallel to the 2PWM at a distance L/2, 
as sketched in Fig. 0. The x-PMTM wave function of momentum K reads: 



whereas (0,0), b= (3,0),c= (0,3), and 



1 

3 



1 =ir,w^ (8) 



(K x - k x ) ■ 6/(2vr) must be odd, which leads to k x = 0,±2?r/3 for K = (2?r/6, 2?r/6). The 
y-oriented PMTM wave function \^ y (K,k y )) is defined in a similar way. The final PMTM 
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Fig. 3. - Upper left: GS projections over PMTM ansatz with k x and k y = 0; Upper right: Relative 
errors AE/E; The bare ansatz behaviors, corrected by a first order or a second order t/U expansion 
of 2PWM are given by dashed lines, dashed-dotted lines and the dotted lines respectively; Lower left: 
Average inter-particle spacings (d); Lower right: Fluctuations Ad of the inter-particle spacings as a 
function of r s . Exact behaviors (thick lines), and ansatz behaviors corrected by a second order t/U 
expansion of 2PWM. K = (2tt/6, 2tt/6). 



ansatz of momentum if is a normalized combination of the x and y-PMTMs, which reads: 

\* X (K,k x ))-\*y(K,ky)) ^ 



\V(K,k x ,ky)) 



^/2-2(V x (K,k 3x ) | %{KMv)) 



! -(\^ x (K,k x )}-\y y (K,k y )}), (9) 



and the constraint k x = k y makes it invariant under x — y permutation. 

In Fig. 2, the three values P X (K, k x ) = \ (^q(K) | $> x (K, k x )) | 2 taken by the projections of 
the exact GS \^>o(K)} over the x-PMTMs \fy x (K,k x )) are given as a function of r s , together 
with the GS projection P(K, k x ,k y ) = \ (^o(K) \ ^(K, k x , k y )) | 2 over the PMTM of momenta 



i k y ) — 

(K,k x — 0, k y = 0). Following the three projections over the x-PMTMs of different wave 
vector k x , one can see how the third particle gets progressively localized in the x-direction as 
r s increases, the rigid three particle triangular molecule corresponding to P X (K, k x ) = 1/3 for 
the three possible k, 
with k 



P(K, k x = 0, k v = 0) w 93% at r s w 100. However, only the PMTMs 
< 50. For those values of r s , the third particle is fully 



k y = contribute when r s 



delocalizcd in the direction parallel to the oriented PMWMs. However, it is likely that the 
PMTM ansatz overestimates the rigidity of the remaining 2PWM when r s becomes smaller. 
This can be partly fixed using a t/U expansion for the 2PWM (as sketched in Fig. 2 left) and 
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keeping the third particle in its delocalized plane wave state with k x (k y ) = 0. 

The improvements coming from this partial t/U expansion of the PMTM ansatz are given 
in Fig. 3, where one can see the behaviors of the bare ansatz, of the ansatz corrected to 
first order and to second order of the t/U expansion of the 2PWM. In the upper figures, 
the GS projections and the relative errors AE(p)/E are shown, E denoting the exact GS 
energy, AE(p) — Ea{p) — E, E^(p) being the ansatz energy at the p th order of the partial 
t/U expansion. Not only the GS description is improved, but lower values of r s can now be 
reached. In the lower figures, the three GS interparticle spacings d rn i„ di nt and d max are given, 
and compared to the corresponding values of the ansatz, after a second order t/U expansion 
of the 2PWM. As underlined by the arrows, both the averages and the variances are now well 
described for r s ss 40. 

However, let us underline that the mesoscopic melting process is not yet achieved at r s ss 40. 
From a study of the weak coupling limit, we have obtained precursor behaviors of the formation 
of the Wigner molecule at smaller r s . For instance, certain GS projections over low energy 
non interacting states which are close in energy to the non interacting GS, but orthogonal 
to the large coupling Wigner molecule begin to decay when r s > 5. Therefore, the PMTM 
ansatz, even improved by a t/U expansion of the 2PWM fails to describe this precursor regime 
(5 < r s < 30) where a floppy 2PWM takes place, but is not rigid enough to be described by a 
simple t/U expansion. 

In summary, we have shown that the quantum melting of a three particle Wigner molecule 
confined on a 2d torus proceeds via an intermediate regime which can be described by the 
simple concept of a partially melted Wigner molecule, built of a delocalized particle and of 
a floppy 2PWM. This is in agreement with the general multi stage picture of mesoscopic 
quantum melting given by other works using 2d parabolic traps. At a mesoscopic scale, this 
gives a simple illustration of the quantum crystal with k = defectons conjectured by Andreev 
and Lifshitz. Notably, one can see that the number of Wigner lattice sites is smaller than the 
total number of charges. This shows that the multi stage melting is not a mesoscopic surface 
effect and suggests that dilute 2d electron gases of intermediate r s could be more complicated 
than usually assumed. 

This work was supported by the EU program "Nanoscale dynamics, coherence and compu- 
tation" and the Hungarian Science Foundation OTKA T025866 and T034832. 
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